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■  Final  Technical  Report  on  FY90  SORP 

i<mo 

Randall  S.  Jacobson 
ONR  Code  1125GG 


My  FY90  SORP,  "Attenuation  of  compressional  waves  in  oceanic  crust",  has  been  fruitful.  The 
proposal  was  to  analyze  a  high  quality  seismic  refraction  data  set  to  obtain  compressional  wave  attenuation 
models  for  oceanic  crust.  The  specific  tasks  and  results  are: 

Complete  the  attenuation  analysis  and  publish  the  results.  Despite  the  late  delivery  of  the  computer  to 
do  the  analysis,  I  was  able  to  finish  the  work  and  write  a  manuscript  (R.  S.  Jacobson  and  B.  T.  R. 

Lewis,  The  First  Direct  Measurements  of  Compressional  Wave  Attenuation  in  the  Uppermost  Oceanic 
Crust).  This  paper  has  been  published  in  the  Journal  of  Geophysical  Research  (v.  95,  Bll,  pp.  17,417- 
17,429)  and  is  attached.  The  basic  conclusions  of  this  work  were: 

(a)  Values  of  Q  (a  measure  of  the  attenuation)  cluster  between  20  and  50,  an  order  of 
magnitude  greater  attenuation  than  estimated  by  forward  modeling  by  others. 

(b)  Variability  of  the  results  strongly  implicate  heterogeneities  within  the  crust  on  scales  of 
lOO’s  of  meters. 

(c)  The  attenuation  estimates  are  considered  normal  to  high,  based  upon  the  use  of  the  sonar 
equation  to  match  signal  to  noise  ratios  as  a  function  of  frequency  of  historical  data. 

(d)  Values  of  Q  around  50  are  sufficient  to  limit  propagation  of  seismic  energy  within  the 
crust  to  be  below  20  Hz  for  typical  seismic  refraction  profiles. 

Additional  work  was  directed  towards  joint  inversion  of  travel  time  and  amplitude  data  for  velocity-depth 
information.  I  have  written  software  to  correct  all  the  seismic  data  for  input  into  the  inversion  program.  The 
development  of  the  software  for  the  inversion,  testing  of  the  inversion  methodology,  and  the  analysis  of  the 
results  has  been  completed.  The  results  were  somewhat  disappointing,  perhaps  due  to  errors  in  location  of  the 
shot  positions  in  three  dimensions.  A  further  test  of  this  methodology  is  expected  in  the  future  by  inversion  of 
another  data  set,  using  electrically  detonated  explosive  shots,  providing  accurate  origin  times  of  die  source. 
Other  progress  relates  to  the  unexpected  task  of  determining  the  rank  of  matrices  for  inversion.  A  completely 
objective  method  has  been  found  to  determine  which  eigenvalues  are  non-zero  using  Chi-square  analysis  and 
Fisher  Tests.  This  combination  of  error  analysis  and  statistical  analysis  has  yet  to  be  applied  to  inversion 
theory  but  is  often  used  for  other  types  of  problems.  I  have  instituted  a  literature  search  for  applications, 
theoretical  developments,  etc.  A  software  package  has  been  developed  to  use  this  methodology,  and  I  am 
currently  using  it  in  inversion  applications. 
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The  First  Direct  Measurements  of  Upper  Oceanic  Crustal 
Compressional  Wave  Attenuation 

R.  S.  Jacobso.v 

Office  of  Nctjc! Research.,  Arlington,  Virginia 


B.  T.  R.  Lewis 

School  of  Oceanography,  University  of  Washington,  Seattle,  Washington 

The  first  Jirecc  measurement  of  compressicnal  wave  attenuation  of  the  uppermost  650  m  of  oceanic  crust  was 
performed  using  data  recorded  by  seafloor  hydrophones  and  large  (56-1 16  leg),  deep,  explosive  sources.  Tne  site 
was  13  lan  east  of  the  southernmost  Juan  de  Fuca  Ridge  cn  crust  0.4  m.y.  old.  Spectral  ratios  were  performed 
between  bottom  refracting  waves  and  direct  water  waves,  adjusted  for  spreading  losses  and  transmission  coefficient 
losses.  Several  tests  of  the  data  were  performed,  demonstrating  that  attenuation  is  linearly  related  to  frequency 
between  15  and  140  Hz,  but  frequency-independent  components  of  attenuation  are  also  evident.  Values  of 
compressional  wave  Q  cluster  between  20  ar.d  50  and  do  not  show  any  systematic  variation  with  depth  over  650 
m.  The  attenuation  results  also  indicate  the  presence  of  heterogeneities  within  the  crust,  as  the  solutions  for  each 
receiver’s  data  set  are  significantly  different.  No  evidence  for  azimuthal  variations  of  attenuation  are  supported 
by  the  data,  although  the  data  do  not  optimally  sample  a  wide  variation  of  azimuths.  Our  attenuation  values  are 
judged  to  be  normal  to  higher  than  expected  for  the  whole  oceanic  crust,  based  upon  comparisons  to  results  from 
synthetic  seismogram  modeling  by  others  and  by  modeling  signal  to  noise  ratios  of  typical  seismic  refraction 
profiles.  The  results  are  consistent  with  recent  laboratory  measurements  at  ultrasonic  frequencies  for  dry  and 
saturated  basalts  at  seafloor  pressures  and  temperatures.  /’ 


Introduction 

Despite  decades  of  investigations  of  the  oceanic  crust  and  mantle 
using  seismic  refraction  techniques,  knowledge  of  the  velocity  and 
attenuation  structure  of  the  uppermost  kilometer  is  the  most  defi¬ 
cient  [ Spudich  and  Or  cast,  ly80a).  Despite  the  fact  that  seismic 
energy  must  propagate  through  this  portion  of  the  crust  to  be 
received  by  either  on-bottom,  midwater  or  near-surface  sensors, 
nature  appears  to  conspire  against  investigations  of  this  most  acces¬ 
sible  portion  of  oceanic  lithosphere.  To  isolate  the  effects  of  the 
uppermost  crust  upon  seismic  wave  propagation  requires  a  con¬ 
comitant  isolation  of  other  interfering  energy  paths.  To  optimize 
this  effort,  both  sources  and  receivers  should  be  near  the  seafloor, 
so  that  interfering  water-borne  energy  will  arrive  at  times  sufficient¬ 
ly  distinct  from  those  paths  traveling  within  the  bottom.  Two 
decades  of  refraction  experiments  using  surface  charges  and  ocean 
bottom  seismometers  have  revealed  little  of  the  detailed  nature  of 
the  uppermost  kilometer  of  crust,  although  some  gross  averages  can 
be  obtained  [Ewing  and  Purdy,  1982],  Recent  multichannel  studies 
with  high  spatial  resolution  [e.g„  Harding  et  al.,  1989;  Vera,  1989] 
have  revealed  details  of  the  uppermost  velocity  structures  and 
seafloor  velocities.  Vertical  or  oblique  seismic  profiling  has  also 
proved  effective  for  determining  upper  crustal  velocities  {e.g., 
Siephenetal,  1 980;  Stephen  and  Harding,  1983;  Little  andStephen, 
1985].  Only  recently  have  experimental  techniques  been  Cried  to 
use  near-bottom  charges  for  seismic  refraction  profiling  (e.g,, 
Jacobson  et  al.,  1981;  Purdy,  1986;  Sauter  et  al,  1986],  with 
success.  There  are  several  disadvantages  using  deep  explosive 
charges,  including  higher  frequency  content,  more  complicated 
logistics,  more  time  required  to  perform  the  experiment,  and  in¬ 
creased  safety  concerns.  The  goal  of  obtaining  more  detailed  seis¬ 
mic  information  of  seismic  layer  2  is  closer  to  being  achieved. 
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particularly  in  concert  with  results  from  drilling  the  oceanic  crust  by 
the  Deep  Sea  Drilling  Project  and  the  Ocean  Drilling  Project 

The  seismic  results  to  date  are  entirely  in  terms  of  compressional 
wave  velocity  structures,  often  showing  steep  gradients  with  depth 
that  are  a  function  of  crustal  age  [Purdy.  1987].  Other  information 
regarding  seismic  structures,  such  as  compressional  wave  attenua¬ 
tion,  shear  velocity,  and  attenuation  structures,  are  either  unavail¬ 
able  for  analysis  in  existing  data  sets  or  have  not  been  investigated. 
Tnese  additional  seismic  structures  can  be  useful  not  only  for  further 
resolving  the  geologic  stratigraphy  but  arc  also  crucial  to  under¬ 
standing  the  full  impact  of  upper  oceanic  crustal  structures  on 
seismic  wave  propagation  into  and  out  of  the  oceanic  lithosphere. 

In  1985,  we  conducted  a  deep  source,  deep  receiver  seismic 
refraction  experiment  on  0.4  m.y.  old  crust  13  km  east  of  the 
southernmost  Juan  de  Fuca  Ridge  (Figure  1).  The  site  was  chosen 
to  coincide  with  a  previous  experiment  that  used  a  deep-towed 
seismic  streamer  and  air  guns  along  a  profile  shot  parallel  to  the 
spreading  direction  [Lewis  and  Jung,  1984],  Analysis  of  those  data 
demonstrated  that  the  frequency  dependence  of  attenuation  was 
essentially  independent  of  observational  range  or  depth  of  penetra¬ 
tion.  Further  analysis  by  Lewis  and  Jung  [1989]  found  that  P  wave 
Q  between  10  and  50  for  the  upper  1  km  at  this  site  were  consistent 
with  the  data.  They  also  suggested  the  possibility  of  velocity 
structures  within  the  bottom  that  create  ’’resonances”  due  to  ducting 
of  seismic  energy  that  are  frequency  dependent  and  modify  any 
attenuation  that  may  occur.  The  experimental  geometry  of  the 
deep-lowed  streamer  data  set  was  suboptimal  to  determine  the  fine 
structure  of  attenuation  parameters  with  depth. 

Experimental  Configuration  and  Design 

To  test  the  hypothesis  for  the  anomalous  attenuation  due  to 
velocity  "resonances"  required  an  unusual  experimental  configura¬ 
tion.  The  goals  were  to  determine  simultaneously  the  comprcs- 
sional  wave  velocity  and  attenuation  structure  with  depth  and  to 
examine  the  azimuthal  dependence  of  attenuation. 

Deep  explosive  charges  were  needed  with  sufficient  energy  and 
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Fig.  I.  A  map  derived  from  Sea  Beam  bathymetry  of  the  experimental  rile  located  13  km  east  of  the  southernmost  Juan  de  Fuca 
Ridge.  The  solid  triangles  indicate  the  locations  of  the  ocean  bottom  seismometers  (1-6)  and  the  seafloor  hydrophones  (7-10);  only 
the  SFH  data  were  used  to  invert  for  attenuation.  The  solid  dots  are  the  shot  locations,  and  the  interconnecting  lines  indicate  the 
shot-receiver  combinations  used  in  the  analysis.  The  ray  paths  are  heavily  biased  to  the  north  and  south,  a  result  of  experimental 
geometry  and  possibly  propagation  effects.  The  contours  are  in  uncorrected  meters. 


bandwidth  to  overcome  the  high  attenuation  expected.  TNT  andC4 
plastic  explosives  were  detonated  by  (SUS)  Sound  Underwater 
Source  [Naval  Ordnance  System  Command,  1973]  MK  59  devices, 
set  to  1829  or  2438  nominal  meters  (6000  or  8000  feet).  Charge 
sizes  ranged  from  58  to  116  kg,  yielding  maximum  energies  at  50-60 
Hz.  Bottom  charges,  designed  to  excite  shear  waves,  consisted  of 
approximately  3  kg  of  explosives  in  a  "pipe  bomb".  These  charges 
were  largely  unsuccessful,  as  the  energy  needed  to  release  the  end 
cap  was  insufficient  to  overcome  ambient  hydrostatic  pressures. 

Receivers  consisted  of  digitally  recording  seafloor  hydrophones 
[Shor,  1979],  with  126  dB  of  dynamic  range,  66  dB  of  resolution, 
and  200  Hz  effective  bandwidth  with  1-ms  sampling.  Due  to  power 
limitations,  two  hydrophone  systems  weTe  deployed  twice;  on  the 
second  deployment,  one  of  the  seafloor  hydrophones  (SFH9)  was 
deployed  approximately  700  m  above  the  bottom,  so  that  uncon¬ 
taminated  source  signatures  could  be  recorded.  In  addition,  six 
triaxial  ocean  bottom  seismometers  were  deployed,  which  recorded 
analog  seismic  signals  on  magnetic  tape,  yielding  a  bandwidth  of 
approximately  100  Hz  and  a  dynamic  range  of  80  dB  after 
decompression  [Johnson  et  ai,  1977].  The  ocean  bottom  seis¬ 
mometer  data  were  not  analyzed  in  this  paper  due  to  the  limited 
bandwidth  of  the  instruments. 


All  of  the  instruments  were  deployed  in  a3-km-diameter  array  in 
an  area  of  fiat  seafloor,  2600  m  deep,  13  km  east  of  the  Juan  de  Fuca 
Ridge  (Figure  1).  Box  and  gravity  cores  recovered  up  to  039  m  of 
fine  dark  brown  silt  over  heavy  grey  clay.  Shipboard  and  bottom 
pingers  were  unable  to  resolve  sediment  thickness,  which  is  thought 
to  be  no  more  than  a  couple  of  meters  thick.  At  the  seismic 
frequencies  of  interest  (  5  150  Hz),  the  sedimentary  cover  should 
have  an  insignificant  impact  upon  experimental  results. 

Navigation  consisted  of  LORAN  C  and  transit  satellites,  both 
generally  accurate  to  100  m  when  available.  Bathymetry  consisted 
of  3.5-  and  12-kHz-wide  beam  echo  sounders  and  Sea  Beam 
coverage  of  the  area  was  provided  by  NOAA. 

Two  separate  sequences  of  shots  were  deployed  in  a  star  pattern 
to  permit  varying  azimuths  and  ranges  to  the  receivers.  Corrections 
to  the  hydrophone  refraction  data  were  reported  by  Poujol  and 
Jacobson  [1988].  All  shots  were  corrected  to  the  entry  point  at  the 
seafloor,  defined  by  the  Sea  Beam  bathymetry  using  the  water  path 
correction  of  Purdy  [1982].  A  record  section  was  plotted  (Figure 
2),  and  travel  times  were  picked  manually.  The  resultant  travel  time 
plot  is  shown  in  Figure  3  with  a  third-order  polynomial  fit,  which 
was  later  used  to  facilitate  the  determination  of  the  velocity  depth 
function  [Poujol  and  Jacobson,  1988]. 
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Attenuation 

To  determine  the  attenuation  parameters  from  seismic  refraction 
data,  one  needs  to  account  for  all  mechanisms  that  can  modify  the 
amplitude  and  frequency  content  of  the  propagating  wave.  The 
received  signal,  A(x,co),  is  affected  by  spreading  losses,  reflection, 
and/cr  transmission  losses,  instrumental  response,  and  attenuation: 

A(*,co)  =  A^O.ca)^*-^  (1) 

L'W 

where  x  is  the  distance,  to  is  the  angular  frequency,  A0(0,co)  is  the 
source  function,  exp[-ct(co)x]  is  the  attenuation  due  to  propagation, 
R(x,cd)  is  the  total  loss  due  to  reflections  and  transmissions,  C{x) 
represents  the  spreading  loss,  and  /(to)  is  the  instrument  transfer 
function.  The  direct  water  wave  path  can  be  expressed  as 

A4r.»)  =  AofO.ta)^^^  (2) 

where  exp[-a»xl  is  the  attenuation  through  the  water  (hereafter  crw 
is  assumed  to  be  zero),  and  subscript  w  represents  the  direct  water 
wave.  Dividing  equation  (1)  by  equation  (2)  and  taking  logarithms, 
we  obtain 

= + MW + (3) 

We  will  assume  that  there  are  r.o  caustics  and  that  all  frequency 
dependence  is  in  a;  that  is,  both  R{x,( o)  and  spreading  losses  are 
independent  of  frequency.  The  term  otr  can  be  rewritten  several 
ways  [e.g.  Johnston  and  Toksdz,  1981],  including 

ct(co)x  =  n  FxQ-\F)V\F)  =  nFT(F)Q-\F)  (4) 

where  T  is  the  travel  time  through  the  attenuating  medium  and  F  is 
the  frequency  in  hertz.  Equation  (4)  allows  for  velocity  dispersion 
and  3  frequency-dependent  quality  factor  Q.  In  general,  although 
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Fig.  3.  Seafloor  hydrophone  travel  time  data  (circles)  used  for  the  inversion 
of  the  velocity  depth  function,  which  in  turn  is  used  to  constrain  the  spreading 
losses.  A  third -order  polynomial  fit  to  the  data  is  also  shown. 

Q  can  also  vary  with  depth,  we  will  assume  for  the  moment  that 
Q{F)  is  constant  throughout  the  bottom.  Given  these  assumptions 
and  definitions,  equation  (3)  can  be  rewritten  as 

(5) 

Considering  only  those  terms  that  are  frequency  dependent  and 
further  assuming  that  T{F)  is  independent  of  frequency,  the  specific 
quality  factor  Q'1  can  be  found  to  be 

Q~\F)  =  [ In f (6a) 

*  v  ;  nFT  [^A(xJr)  J 


Distance,  km 

Fig.  2.  True  amplitude,  unreduced,  record  section  of  the  seafloor  hydrophone  dau  with  the  predicted  travel  lime  curve,  based  upon 
a  third  order  polynomial  The  traces  were  terminated  at  the  water  wave  arrival  and  equalized  for  shot  size  variations  by  where 
W  is  the  weight  in  pounds.  No  clear  evidence  of  converted  shear  waves  is  present  in  the  data.  The  chaotic  nature  of  the  imenrace 
coherence  is  due  to  the  widely  varying  shot  depths,  which  affects  both  the  frequency  content  of  the  source  function  and  the  amplitudes, 
due  to  spreading  losses. 
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If  the  logarithm  of  the  spectral  ratio  is  assumed  or  found  to  be 
linearly  dependent  upon  frequency  and  if  one  assumes  that  the 
frequency-independent  terms  of  equation  (5)  do  not  affect  Q  in 
any  way,  then 


where  Q  l  is  now  independent  of  frequency.  The  assumptions  of 
frequency-independent  Q ‘l,  T,  and  a  linear  dependence  of  the 
logarithm  of  the  spectral  ratio  with  frequency  violate  the  principles 
of  causality  [cf.  Kjartensson,  1979  and  Jacobson,  1937],  Equation 
(6b),  although  physically  unrealizable,  is  often  used  as  the  spectral 
ratio  method  for  seismic  refraction  profiles  [e.g.,  Jacobson  et  a!., 
1981,  1984],  Equation  (6j)  is  the  more  appropriate  relationship  to 
investigate  attenuation,  despite  the  fact  that  Jacobson  [1987] 
demonstrated  that  velocity  dispersion,  V(F),  or  equivalently  T(F), 
as  well  as  frequency-dependent  Q  can  be  measured  using  seismic 
refraction  methodologies.  Frequency  dependence  of  the  attenuation 
coefficient  can  be  tested  simply,  using  the  form 

a(F)  =  xF”TQ~x  (7) 

lumping  the  frequency  dependence  of  T  and  Q' 1  into  the  power  of 
frequency  dependence,  m.  Attenuation  due  to  scattering  can  also 
exhibit  frequency-ir.dependentcomponents,  so  thatabetter  testmay 
be 

a(F)  =  nFmTQ;l  +  B  (8) 

where  B  is  a  term  independent  of  frequency.  Jacobson  [1987] 
expanded  this  definition  to 

a(/0  =  JtfT^;1f"  +  y  +  yJ-A3  (9) 

to  search  for  a  more  complete,  empirical,  description  of  the  effects 
of  scattering  attenuation.  To  do  so  requires  the  calculation  of  the 
R(x)  and  G(z)  terms,  which  can  be  misconstrued  as  frequency-inde¬ 
pendent  components  of  attenuation. 

The  spreading  loss  term,  G(x),  for  a  flat  earth  can  be  derived  from 
the  equation  for  a  spherically  symmetric  earth  [Aki  and  Richards, 
1980,  p.  120]: 

G(x)  =  y|cos(9,)  cos(0,)  ^  (10) 


where  p  is  the  slowness,  9  refers  to  the  angle  of  the  ray  with  respect 
to  the  vertical  and  subscripts  r  and  s  refer  to  the  receiver  and  source, 
respectively.  The  term  dx/dp  equal  to  (d*l/dx?)  1  can  be  obtained 
from  the  travel  time  curve  or  by  manipulation  of  the  velocity  depth 
function  to  yield  xip).  In  a  layered  media,  x  must  be  replaced  by  Lx 
and  dxjop  by  Zbxldp  over  the  layers  the  wave  propagates. 

The  reflection  and  transmission  coefficients  can  be  calculated 
from  equations  given  by  Aki  and  Richards  [1980],  if  one  knows  or 
assumes  the  compressional  and  shear  wave  velocities  and  densities 
throughout  the  entire  structure  of  interest.  Further  assumptions 
include  planar  waves,  smooth  boundaries,  and  no  velocity  gradients, 
all  of  which  can  modify  the  reflection  and  transmission  coefficients. 
The  appropriate  coefficients  can  only  be  predicted  with  reasonable 
accuracy  ax  the  seafloor  interface.  Since  the  details  of  the  subbottom 
structure  are  unknown,  we  are  forced  to  incorporate  subbottom 
losses  due  to  transmissions  and  reflections  into  the  attenuation 
coefficients. 

Equation  (9)  can  be  used  to  examine  the  data  for  frequency-de¬ 
pendent  quality  factors,  Q 1 ,  as  well  as  frequency-independent  terms 
that  indicate  a  “DC"  shift  in  attenuation.  Jacobson  [1987]  has 
interpreted  these  terms  to  indicate  either  a  frequency  dependence  of 


Q  below  the  lowest  frequency  measured  or  offsets  of  attenuation 
due  to  interbed  multiples  or  to  scattering  from  heterogeneities. 

Provided  the  data  indicate  no  frequency  dependence  of  Q  within 
the  measured  bandwidth  and  relaxing  the  assumption  of  a  depth 
independent  Q,  we  can  solve  for  Q(z).  To  do  so  requires  knowledge 
of  the  ray  path  geometry  for  each  shot,  easily  obtained  by  ray  tracing, 
knowing  the  velocity  depth  function.  This  method  allows  a  direct 
mapping  of  total  travel  time  T  or  total  subbottom  distance*  into  an 
equivalent  bottoming  depth  z  where  the  ray  turns  horizontal  [cf. 
Dorman  and  Jacobson,  1981].  Keeping  the  frequency-independent 
terms,  the  more  general  form  of  equation  (5)  becomes 

=  -rJT2_’ (z)  +  hlR{x)  +  *{  CW  ]  0 1 } 

Assuming  an  attenuation  model  consisting  of  aseries  of  layers,  each 
with  a  constant  Q'\  and  retaining  the  discretized  frequency,  F,, 
equation  (1 1)  can  be  rewritten  as 


02i 


where  subscript  j  denotes  a  particular  shot-receiver  pair  and  sub¬ 
script  k  denotes  the  layer  number.  Tjt  represents  the  travel  time  of 
the  ;th  ray  in  layer  k  similarly,  Q a’1  is  the  contribution  of  attenuation 
to  the  spectral  ratio  for  layer  k  equal  to 

Q'^Q-ol  +  ~  (13) 

fi 


where  B  describes  the  frequency-independent  components  of  at¬ 
tenuation  a  analogous  to  equation  (8). 

Once  the  effects  of  attenuation  are  isolated  from  all  other 
parameters,  the  azimuthal  dependence  of  attenuation  can  be  deter 
mined.  An  empirical  test  of  azimuthal  dependence  of  attenuation 
can  be  performed  via 

g"1  (;,(?)  =  gri(:)X  a,cos(<p,)  +  b,sin(<?,)  (14) 


Only  even  powers  of  a  and  b  should  exist,  as  the  presence  of  odd 
powers  would  violate  reciprocity. 

Velocity  Results 

To  analyze  properly  the  data  for  attenuation  as  a  function  of  depth 
it  is  necessary  to  determine  the  velocity  depth  function.  Poujol  and 
Jacobson[  1988]  reported  on  the  velocity  inversion  for  the  seafloor 
hydrophone  data  only,  using  the  x-C,  inversion  of  Dorman  and 
Jacobson  [1981].  Their  results  indicate  a  seafloor  velocity  of  2.7 
km  s'1,  wiLh  high  velocity  gradients,  ranging  from  4.6  s'1  at  the 
surface  to  4.1  s'1  at  679  m  depth,  the  maximum  extent  of  the  data. 
The  maximum  velocity  observed  was  5.6  km  s' 1 .  The  velocity  depth 
profile  (Figure  4)  is  extremely  smooth,  showing  no  sharp  velocity 
changes  usually  associated  with  layer  2A-B  stratigraphy.  This  is 
probably  not  real  but  is  simply  the  result  of  oversmoothing  of  the 
travel  time  curve  used  to  invert  the  velocity  depth  function. 

Poujol  and  Jacobson  [  1 988]  point  out,  however,  that  considerable 
scatter  in  their  waveforms  may  be  due  to  local  crustal 
heterogeneities.  In  particular,  one  seafloor  hydrophone  (SFH8) 
recorded  an  unusual  high-frequency,  low-amplitude  precursor  on 
several  of  the  records,  due  to  what  might  be  a  local  high-velocity 
block  of  crustal  material.  Further  evidence  of  crustal 
heterogeneities  can  also  be  found  in  the  ocean  bottom  seismometer 
data.  Although  a  similar  inversion  technique  for  velocity  structure 
using  all  available  data  [Jung,  1988]  revealed  similar  structures  to 
those  as  reported  by  Poujol  and  Jacobson  [1988],  two  of  the  ocean 
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Fig.  4.  The  velocity  depth  solution  with  95  7o  error  bounds  and  the  velocity 
gradient  solution  obtained  by  inverting  the  travel  time  data. 


bottom  seismometers  (OBS2  and  OBS6)  recorded  converted  shear 
waves.  Assuming  a  Vp/V,  ratio  of  1.8  at  the  seafloor,  a  shear  wave 
velocity  of  1 .5  km  s" 1  would  permit  conversion  to  shear  wave  energy 
at  the  seafloor.  If  the  shear  velocity  dips  below  the  overlying  water 
velocity,  conversion  to  shear  energy  becomes  increasingly  ineffi¬ 
cient  If  the  shear  (and  compressional)  velocities  at  the  seafloor  vary 
by  0.1  km  s'1  due  to  crustal  inhomogeneities,  significant  variations 
of  shear  wave  energies  will  result  No  clear  converted  shear  wave 
energy  was  observed  on  the  hydrophone  data  (Figure  2). 

For  the  attenuation  analysis,  only  the  seafloor  hydrophone  data 
will  be  used,  as  the  ocean  bottom  seismometer  data  were  too 
band-limited.  To  be  internally  consistent,  the  velocity  structure 
determined  by  Poujol  and  Jacobson  [1988J  will  be  used  in  the 
subsequent  analysis.  The  differences  between  the  velocity  struc¬ 
tures  determined  by  Jmg  [198  8]  and  Poujol  and  Jacobson  [1988] 
are  minor;  the  effects  on  the  attenuation  results  using  either  velocity 
structure  are  insignificant. 

Attenuation  Analysis 

To  determine  the  attenuadon  structure  versus  depth,  one  must  first 
calculate  the  spectral  ratio  (equadon  (6)).  Both  the  refracted  energy 
and  the  direct  waterborne  arrivals  were  box-car  windowed  into 
0.064-s  segments,  except  the  most  close-in  arrival,  which  was 
windowed  into  a  0.032-s  time  segment  No  other  smoothing  opera¬ 
tions  were  performed  on  the  data.  Only  the  first  refracted  arrivals 
were  used;  no  other  refracted  phases  could  be  identified.  Fast 


Fourier  transforms  of  the  windowed  data  were  performed,  the 
spectral  rados  determined,  and  converted  into  logarithms.  Typical 
spectra  (FigureS)  displayed  a  flattening  of  the  spectral  levels  beyond 
150  Hz,  indicating  digital  noise  due  to  limited  resolution.  Each 
spectral  rado  was  examined  for  sufficient  signal  to  noise  rados,  and 
only  data  with  adequate  signals  were  retained  for  further  analysis. 


Fig.  5.  (Top)  A  seismic  trace  recorded  by  SFH9  at  a  reduced  distance  of 
2.154  km,  beginning  at  the  first  refracted  arrival  time.  The  direct  water 
wave,  which  dominates  the  trace,  was  windowed  with  a  0.064-s  box-car  for 
spectral  analysis.  The  vertical  axis  is  in  digital  counts.  (Middle)  The  same 
aeismic  trace  as  above,  except  expanded  to  show  only  the  refracted  arrivaL 
There  is  a  small  DC  offset  in  the  trace  due  to  the  gain-ranging  amplifiers. 
The  first  64  ms  of  data  were  selected  for  spectral  analysis.  (Bottom)  The 
ipectral  plot  of  the  direct  water  wave  (large  dashes),  the  refracted  wave 
(small  dashej)  and  the  corrected  ipectral  ratio  (solid).  The  ordinate  is  in 
decibels  relative  to  digital  counts  for  the  direct  and  water  wavet.  The 
ipectral  ratio  data,  adjusted  for  spreading  losses  and  transmission  loasei,  are 
in  decibels.  The  circles  represent  the  data  actually  used  for  the  inversion  of 
attenuation. 
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Most  of  the  data  *e:e  coofincd  to  the  bandwidth  15-94  Hz,  although 
some  dau  extended  to  140  Hz.  To  examine  the  existence  of  fre¬ 
quency-independent  components  of  attenuation,  the  raw  spectral 
ratio  data  w  ere  adjusted  for  spreading  losses  using  equation  (10)  and 
for  transmission  losses  across  the  water-basement  interface. 

Tests  for  Frecuenc-  Dependence  of  Attenuation 

The  first  examination  of  the  data  was  to  look  for  frequency 
dependence  cf  the  attenuation  coefficient  a{F)  by  using  equation 
(9).  Various  permutations  of  the  coefficients  were  examined,  as  the 
coefficients  are  not  orthogonal.  F  tests  were  used  to  determine  the 
significance  of  variance  reductions  with  the  addition  of  each  new 
coefficient.  Results  of  the  tests  indicate  that  the  combined  data  set 
(ALL)  should  be  broken  down  into  individual  data  sets  appropriate 
to  each  of  the  four  receivers,  as  discussed  below.  In  all  five  cases, 
the  attenuation  coefficient  was  found  to  have  a  linear  dependence 
of  frequency  (r.=0  of  equation  (9))  within  the  bandwidth  of  available 
data  (15.625-140.625  Hz).  Each  data  set  did  exhibit  some  dme- 
ar.d/or  frequency-independent  components  of  artenuadon  (e.g..  At, 
A;  and  A3  of  equation  (9)).  In  three  cases  (SFH8,  SFH10,  and  ALL) 
some  frequency-independent  component  of  attenuation  was  ex¬ 
hibited  (Tafcie  1).  These  coefficients  can  be  interpreted  as  a  non¬ 
zero  intercept  of  a  linear  fit  to  the  adjusted  spectral  rauo  data  for 
each  shot.  By  calculating  a  slope  and  intercept  for  each  shot,  a 
qualitative  fit  of  the  models  in  Table  1  to  the  data  can  be  seen  in 
Figure  6.  In  all  cases,  the  predicted  slopes  increase  with  bottoming 
depth  (equivalent  to  subbottom  travel  dme  or  range).  For  SFH10 
and  ALL,  the  intercepts  are  largely  negative,  increasing  to  slightly 
positive  values  below  550  m.  Negative  intercepts  for  this  model 
(equation  (9))  imply  an  excess  attenuation  at  zero  frequency,  where 
attenuation  hypothetically  should  be  zero.  Anomalous  sources  of 
attenuation  (e.g.,  scattered  or  apparent  attenuation)  and/or  incorrect 
adjustments  (spreading  losses,  transmission  losses)  are  likely  ex¬ 
planations  for  the  Intercept  offsets. 

SFH7  and  SFH9  models  exhibit  a  component  of  spectral  slope 
other  than  zero  at  the  seafloor  (Figures  6a  and  6c),  resulting  from 
the  time-independent  component  of  attenuation  A2.  Since  equation 
(9)  is  simply  an  empirical  description  of  attenuation,  it  is  difficult  to 
hypothesize  a  source  for  this  component  of  a. 

All  of  the  models  have  values  of  Qo  between  23  and  67,  much 
lower  than  those  values  determinedby  synthetic  seismogT  am  model¬ 
ing  (e.g.,  Spudich  and  Orcutt,  19806],  There  is  little  evidence  of 
range  (or  depth)  independent  variations  of  the  spectral  slope,  since 
most  of  the  data  clusters  between  -0.4  and  -0.8  dB  Hz'1  (Figure  6). 
Similarly,  most  of  the  spectral  intercepts  of  the  data  cluster  around 
zero,  suggesting  the  adjustments  made  to  the  data  are  largely  correct. 
This  latter  observation  can  be  construed  as  evidence  for  the  lack  of 


caustics  in  the  travel  time  data,  which  would  otherwise  produce  a 
pronounced  excess  of  intercept  values  ax  range  or  depth. 

Statistical  tests  (e.g.,  Bartlett’s  test  for  homogeneity  of  variance 
[Ostle  and  Mensing,  1975])  indicate  a  significant  difference  in 
variances  between  the  individual  receiver  data  sets,  implying  the 
differences  in  these  models  are  real.  Additionally,  the  four  in- 
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Fig.  6a.  Plots  of  the  fit  of  the  solutions  using  equation  (9)  and  Table  1  to  the 
various  data  sets,  after  equation  (9)  was  converted  to  decibels  by  multiplica¬ 
tion  of  8.686.  (Top)  The  fil  of  the  spectral  slope  component  of  equation  (9) 
as  a  function  of  bottoming  depth  to  a  regression  of  each  shot’s  spectral  ratio 
data  of  SFH7  by  a  slope  and  intercept.  The  circles  representthebest  fil  slope, 
and  the  error  bounds  correspond  to  a  ±la  range.  The  inversion  for  attenua¬ 
tion  uses  the  data  at  each  frequency  for  each  shot  and  does  not  attempt  to  fit 
the  slope-intercept  combination.  (Bonom)  The  fit  of  the  solution  using 
equation  (9)  to  the  spectral  intercept  of  SFH7  data. 


TA3LE  1.  Model  Solution  of  SFH  Data  Sets  to  Equation  (9),  Before  Conversion  to  Decibels  (Multiplicative  Factor  of  8.686). 


SFH7 

SFH8 

SFH9 

SFH10 

ALL 

Qo' 

0.0215(0.0052) 

0.0307  (0.0067) 

0.0150(0.0174) 

0.0444  (0.0052) 

0.0423(0.0025) 

At 

— 

— 

— 

-1.4789(0.4094) 

-1.9004(0.2244) 

Az 

0  COOS  (0.C022) 

— 

0.0147(0.0091) 

— 

— 

M 

— 

-3.2075(0.4418) 

— 

-2.6350(0.05617) 

-3.1754(0.02930) 

Number  of  tlau 

55 

18 

23 

37 

133 

O 

6.61  dB 

11.15  dB 

5.41  dB 

7.87  dB 

$.73  dB 

The  values  in  parentheses  ire  the  one  standard  deviation  values  of  the  solution  parameters.  The  Urge  dashes  indicate  that  these  parameters  were  not 
statistically  different  from  zero  and  therfore  were  not  required  to  fit  the  data. 
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6b.  Same  u  Figure  6a,  except  using  only  the  data  from  SFHS.  Fig.  6c.  Same  is  Figure  6a.  except  using  only  the  data  from  SFH9. 
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dividual  data  sots  have  a  combined  residual  variance  significantly 
lower  than  the  medal  firing  all  the  data  simultaneously.  It  appears 
that  within  the  3 -km  array,  differences  of  attenuation  of  the  oceanic 
crust  exist  on  this  spatial  scale.  Similar  results  for  spatial 
heterogeneities  of  velocity  within  the  crus:  have  been  determined 
by  Stephen  [19SS]  for  DSDP  site  504B  and  by  Swift  and  Stephen 
[19S9]  for  old  Atlantic  crust. 

Attenuation  Versus  Depth  Results 

Or.e  of  the  conclusions  of  fitting  the  data  to  equation  (9)  was  that 
a  varies  linearly  with  frequency,  so  equations  (12)  and  (13)  can  be 
used  to  examine  the  detailed  aspects  of  variations  of  attenuation  with 
depth.  One  set  of  models  for  the  five  data  sets  consisted  of  a  single 
layer  of  attenuation.  The  results  were  quite  similar  to  those  inTable 
1  and  Figure  6  and  thus  are  r.ot  shown.  Or.e  difference,  however,  is 
that  equations  (12)  and  (13)  force  the  spectral  slope  and  the  special 
intercept  components  to  be  zero  at  the  seafloor.  Only  the  model 
fating  SFHS  data  required  a  nonzero  value  of  the  spectral  intercept 
at  depth.  Values  of  00  of  all  the  models  varied  between  20  and  30, 
clustering  around  the  lower  value.  Again,  there  was  a  significant 
reduction  of  residuals  using  the  individual  data  sets  relative  to  data 
set  ALL,  and  the  variances  between  the  four  individual  data  sets 
were  again  inhomogeneous.  This  emphasizes  the  interpretation  of 
attenuation  heterogeneities  within  the  hydrophone  array. 

The  next  set  of  models  were  determined  by  maximizing  the  depth 
resolution  of  attenuation.  Up  to  14  layers,  matching  the  layering  of 
the  velocity  depth  function,  were  allowed  to  be  fit  with  a  Q  1  and 
the  B'  parameter  (equation  (13)).  Some  of  the  layers  contained 
multiple  rays  turning  horizontally  within  the  layer.  The  data  consist 
of  the  corrected  special  ratio  values  as  a  function  of  frequency,  as 
shown  in  Figure  5.  A  singular  value  decomposition  routine  [Wig¬ 
gins,  1972]  was  used  to  search  for  the  maximum  number  of  eigen¬ 
values  and  eigenvectors  permissible  with  positive  values  of  Q  .  F 
tests  were  used  to  determine  whether  the  reduction  of  variance  with 
increasing  number  of  eigenvalues  was  justified.  The  solutions  to 
the  five  data  sets  are  shown  in  Figure  7  and  chc  fit  to  the  data  in 
Figure  8.  The  data  and  associated  errors  shown  in  Figure  8  are  the 
slopes  and  intercepts  found  by  linear  regression  for  each  shot’s 
special  ratio  values  and  are  used  solely  for  qualitative  examinations 
of  the  model  Fit  to  the  data.  All  of  the  solutions,  except  that  for 
SFH9,  reauired  a  nonzero  intercept  term,  typically  between  15  and 
-15  dB  s"  .  The  values  of  Q0  vary  between  4  and  275  but  cluster 
between  20  and  50.  F  tests  again  indicate  that  the  individual  models 
are  better  than  the  model  fitting  the  ALL  data  set.  Tests  for 
homogeneity  of  variance  indicates  the  variances  are  now 
homogeneous  between  the  various  receiver  models,  suggesting  that 
portions  of  the  individual  solutions  may  have  some  degree  of 
commonality.  The  residual  standard  deviations  for  models  Firing 
data  sets  SFH7,  8,  9,  10,  and  ALL  are  5.63,  7.07,  5.27,  4.71,  and 
7.77  dB,  respectively.  These  residuals  are  much  lower  than  those 
determined  by  fitting  the  data  to  equation  (9),  not  surprising  con¬ 
sidering  the  increased  number  of  parameters  used  to  fit  the  data. 
These  models,  shown  in  Figures  7  and  their  Fit  to  the  data  in  Figure 
8,  arc  our  preferred  solutions  to  the  data. 

Finally,  the  residuals  from  the  solutions  in  Figure  7  were  ex¬ 
amined  for  azimuthal  variations  usmg  equation  (14).  Only  SFH7 
residuals  had  a  significant  reduction  of  variance  using  a  19  term  for 
the  intercept.  This  result  is  physically  unreasonable,  as  it  implies 
reciprocity  does  not  hold.  It  can  be  interpreted,  however,  as  indicat¬ 
ing  some  degree  of  heterogeneity,  since  the  data  set  (as  well  as  all 
of  the  others)  docs  not  truly  sample  all  azimuths  and  there  are  no 
reciprocal  paths. 


Discussion 

Perhaps  the  most  intriguing  result  of  this  study  was  the  consis¬ 
tently  low  values  of  Q ,  regardless  of  which  model  was  used  to  fit 
the  data.  The  values  of  the  quality  far .or  cluster  tightly  between  20 
and  50,  much  lower  than  those  determined  by  Jacobson  eta!.  [1981, 
1984]  for  sediments  but  similar  to  those  inferred  by  Lewis  and  Jung 
[  1 989]  for  oceanic  crust  at  our  site.  V alues  of  Q  obtained  by  forward 
modeling  of  synthetic  seismograms  [cf.  Spudich  and  Orcutt,  1980i>; 
Collins  et  a!.,  1989]  suggest  values  of  Q  closer  to  500,  an  order  of 
magnitude  less  attenuation.  Vera  [1989]  examined  multichannel 
data  from  the  East  Pacific  Rise  and  found  a  Q  of  about  80  to  match 
amplitude  variations  as  a  function  of  offset  Spudich  and  Orcuxt 
[19S0o]  admit  that  their  conslaints  for  Qp  based  upon  synthetic 
seismograms  are  poor,  likely  the  result  of  the  lade-off  of  amplitudes 
between  the  velocity  depth  function  (which  alters  the  spreading  loss) 
and  attenuation.  Further,  the  narrow  band  source  functions  typically 
used  in  synthetic  seismogram  modeling  to  increase  computational 
speed  and  to  mimic  the  resulting  energy  special  content  of  the 
seismic  refraction  profiles  limits  the  frequency  resolution  necessary 
to  consiain  attenuation. 

Wepfer  and  Christensen  [1987,  this  issue]  reponed  the  first 
laboratory  measurements  of  Q  for  dry  and  water-saturated  oceanic 
basalts  under  pressure  and  temperature  conditions  appropriate  for 
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Fig.  7a.  Solutions  of  the  attenuation  varying  as  i  function  of  depth  using 
equations  (12)  and  (13),  after  conversion  to  deribcla,  and  maximizing  the 
depth  resolution  for  SFH7.  (Top)  Plot  of  Q0 ’’  is  a  function  of  depth  (solid) 
and  ±lo  error  bounds  (dashed).  (Bottom)  Plot  of  B',  the  frequency-inde¬ 
pendent  component  of  attenuation,  as  a  function  of  depth  (solid)  and  ±lo 
error  bounds  (dashed). 


17.426 


Jacobson  and  Lewis:  Oujanic  Crustal  Q 


the  oceanic  crust.  Their  results  indicate  Q  between  8  and  100  at 
ultrasonic  frequencies,  depending  upon  the  state  cf  alteration  and 
porosity. 

There  are  no  direct  determinadons  of  values  of  compressional 
wave  Q  for  oceanic  crust  in  the  literature  that  we  are  aware  of.  B  ased 
upon  the  synthetic  seismogram  results  of  others,  we  can  ask  whether 
our  results  are  anomalous  or  typical  for  the  whole  oceanic  crust  To 
answer  this  question,  we  will  assume  a  constant  value  of  Q  for  the 
entire  oceanic  crust  and  apply  equation  (1)  to  a  typical  seismic 
refraction  experiment  using  near-surface  sources  and  receivers.  We 
will  examine  the  signal  to  noise  ratio  for  refracted  arrivals  turning 
horizontally  within  layer  3  at  a  range  where  caustics  and  retrograde 
branches  are  unimportant.  Three  shot  sizes  (13.6, 29. 1.  and  58*2  kg 
or  30,  64,  and  128  pounds,  respectively)  will  be  modeled  as  detonat¬ 
ing  at  quarter  wavelength  depths  [ Shor ,  1963]  as  will  receiver 
depths,  to  maximize  the  constructive  interference  of  sea  surface 
reflections.  The  crustal  velocity  profile  used  will  be  that  of  Fanfare 
4,  as  reported  by  Spudich  and  Orcuit  [19806],  modified  to  eliminate 
the  thin  sedimentary  cover.  Spreading  losses  and  transmission 
losses  were  calculated  for  a  range  of  32.48  km  and  a  slowness  of 
0.14613  s  km'1.  Ambient  noise  levels  were  extracted  from  Urick 
[1986].  The  results,  assuming  a  Q0  of  50  are  shown  in  Figure  9a. 
No  energy  beyond  8  Hz  can  be  received  due  to  the  ambient  noise 
levels.  Maximum  signal  to  noise  ratios  (SNR)  vary  between  1 1  and 
22  dB,  depending  upon  shot  size  and  frequency.  These  SNR  values 
are  consistent  with  the  practices  used  in  two-ship  or  scmobuoy 
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Fig.  8 a.  Fit  of  the  ittcr union  versus  depth  solution  shown  in  Figure  7  to  the 
spectral  data  of  SFH7.  Sec  Figure  6a  for  more  information.  The  predicted 
curves  of  spectral  slope  and  attenuation  are  only  approximately  correct 
between  data  points  and  interpolation  should  be  avoided. 


experiments  where  the  dynamic  range  of  the  recording  xysters  L. 
typically  limited  to  about  20  dB. 

Figure  9 b  shows  the  effect  of  changing  Q  while  keeping  the  shot 
weight  constant  at  29.1  kg  (64  pounds).  Q  of  20  for  the  whole  crust 
would  prohibit  any  recovery  of  refracted  energy.  Q  of  450  would 
exhibit  some  refracted  energy  beyond  20  Hz.  Little  if  any  energy 
was  ever  seen  at  these  ranges  (G.  G,  Shor,  Jr.,  personal  communica¬ 
tion,  1989),  although  some  air  gun  profiles  have  likely  exhibited 
higher  frequency  refracted  arrivals.  A  Q  of  450  would  also  yield 
SNR  values  close  to  35  dB  which  would  be  within  the  range  of 
recorded  levels  using  multiple  recording  gains. 

Based  upon  the  results  of  Figure  9,  it  is  probably  safe  to  say  that 
our  attenuation  results  are  close  to  being  typical  for  the  whole 
oceanic  crust,  although  the  lower  crust  may  have  lower  values  of 
attenuation.  Further  experiments  are  required  to  better  place  our 
results  in  an  overall  picture  of  crusfal  attenuation  values. 

Our  attenuation  results  (Figure  7)  do  not  show  any  systematic 
variation  with  depth,  a  disappointing  result  for  those  interested  in 
correlating  attenuation  with  geologic  stratigraphy.  The  data  set  used 
lies  entirely  within  layer  2,  the  extrusive  basalt  complex.  Neither 
the  attenuation  nor  velocity  profiles  show  any  sublayering  within 
the  first  650  m  of  oceanic  crust  at  this  site. 

The  solutions  of  attenuation  do  indicate  the  possibility  of  a  large 
excess  of  frequency-independent  attenuation  within  the  first  100  m 
(Figures  lb  and  Id).  This  may  be  a  real  result,  due  to  scattered 
attenuation  caused  by  the  rough  seafloor  topography  or  shallow 
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Fig.  86.  Same  u  Figure  8a,  except  uiirg  only  the  data  from  SFH8. 
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g.  Sc.  Same  as  Figure  Su,  except  using  only  the  data  from  SFH9. 
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Fig.  &i.  Same  as  Figure  &j,  except  using  only  too  data  from  SFH10.  In  this 
case,  all  eigenvalues  and  eigenvectors  (see  text  for  details)  were  retained,  so 
that  the  fit  to  the  slope  and  intercepts  are  exact. 


laycring/heterogeneiries.  This  interpretation  is  reinforced  by  no  ting 
that  both  S  FH8  and  SFH10  data  sets  exhibit  this  anomalous  attenua¬ 
tion,  and  both  receivers  are  located  close  to  one  another  (Figure  1). 
Poujol  arulJacobson{\98S ]  reported  anomalous  travel  time  arrivals 
for  SFH8,  suggesting  crustal  heterogeneities  in  the  immediate 
Vicinity  of  the  receiver. 

An  alternative  explanation  for  the  anomalous  frequency-inde¬ 
pendent  attenuation  within  the  first  100  m  may  be  due  to  inap¬ 
propriate  spreading  loss  calculations,  based  upon  the  velocity  depth 
function.  The  method  of  Dorman  and  Jacobson  [1981] 
produces  overly  smooth  velocity  depth  profiles,  due  to  the  averaging 
of  the  travel  time  curve  by  a  low-order  polynomial.  The  spreading 
loss  term  (equation  (10))  can  be  simplified  for  a  constant  velocity 
gradient  layer.  Since 

*  =  ^  cosO  =  (l-pM)  (15; 

?  dz 


where  is  the  velocity  at  the  topandp"1  is  the  velocity  at  the  bottom 
of  the  layer.  Solving  for  the  partial  derivative  of  x  with  respect  to 
slowness  yields 


dx 

dp 
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so  that  the  spreading  loss  term  becomes 
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Fig.  9j.  Plot  of  the  expected  received  seismic  pressure  in  dB  re  1  |aPa  Hz"1^ 
of  a  layer  3  arrival  and  ambient  noise  levels  as  a  function  of  frequency  for  a 
typical  two-ship  seismic  refraction  experiment.  The  range  is  32 1cm,  using 
the  velocity  depth  function  cf  Fanfare  4  [Spudich  and  Orcult,  19806],  The 
received  energy  levels  differ  according  to  shot  size:  13.6  kg  (30  pounds, 
dashed),  29.1  kg  (64  pounds,  solid)  and  53.2  kg  (123  pounds,  dotted)  of  TNT. 
All  of  the  curves  assume  a  Q  of  50  for  the  entire  oceanic  crust.  See  text  for 
more  details. 


Fig.  96.  Same  as  Figure  9a,  except  the  charge  size  was  fixed  at  29.1  kg  (64 
pounds)  and  attenuation  was  varied:  Q  equals  20  (dotted),  50  (dashed)  and 
450  (solid);  the  latter  value  is  normally  assumed  for  oceanic  crust  In  all 
cases,  little  energy  can  be  recorded  above  20  Hz,  in  accord  with  observations 
ever  the  past  three  decades  [Shor,  1963].  While  values  for  Q  of  20  are 
unlikely  for  the  whole  oceanic  crust,  our  results  of  Q  between  20  and  50  are 
largely  compatible  with  values  expected  for  the  oceanic  crust. 


To  increase  the  spreading  less  G(x )  at  shallow  depths,  one  must 
decrease  the  velocity  gradients.  to  account  for  the  excess 

attenuation.  It  is  possible,  therefore,  to  have  a  more  uniform-like 
shallow  layer,  assuming  that  there  is  only  intrinsic  attenuation.  In 
the  same  vein,  the  deficit  cf  attenuation  at  depth  (Figures  8a  and  8 d) 
implies  that  we  have  underestimated  the  velocity  gradients. 

Our  results  also  demonstrate  circumstantial  evidence  for  crustal 
heterogeneities  within  the  3 -km -wide  array.  Every  attempt  to  model 
the  data  statistically  requires  a  separate  solution  for  each  receiver. 
Some  of  the  models  also  required  a  nonphysical  10  dependence  of 
attenuation.  Anomalous  travel  time  residuals  and  arrivals  have  also 
been  noted  by  Poujol  oral  Jacobson  [1988].  Tomographic  aversion 
of  attenuation  would  be  the  logical  next  step,  but  the  pauc  /  of  data 
and  the  few  overlapping  paths  would  render  this  exercise  futile. 

Examination  of  the  data  for  azimuthal  variations,  possibly  due  to 
aligned  vertical  cracks,  showed  no  significant  correlation.  It  must 
be  pointed  out,  however,  that  our  data  set  is  heavily  biased  towards 
ray  paths  in  a  north-south  orientation,  due  to  the  position  or  the 
receivers  relative  to  the  shots  (Figure  1).  Refracted  arrivals  from 
the  west  and  east  were  observed,  but  the  shot  points  were  located 
outside  the  array  and  thus  too  poorly  located  to  be  useful  [Poujol 
ard  Jacobson,  1983],  A  qualitative  comparison  of  an  arrival  at 
SFH3  from  the  east  and  from  the  south  at  approximately  the  same 
range  (Figure  10)  show  a  clear  amplitude  difference.  Thus 
azimuthal  variation  of  attenuation  within  the  shallow  crust  cannot 
be  ruled  out  entirely. 

Finally,  Lewis  and  Jung  [1989]  examined  a  midwatcr  towed  array 
data  set  using  surface  charges  in  the  same  area  as  our  experiment 
Their  forward  modeling  of  attenuation  indicates  that  a  minimum  of 
attenuation  exists  at  10  Hz  due  to  the  presence  of  high-velocity 
gradients  in  the  uppermost  crust,  which  acts  to  duct  high-frequency 
energy  within  the  high  gradient  zone.  Tney  modeled  the  attenuation 
structure  with  two  layers:  a  1  km  thick  upper  layer  with  a  Q  of  50 
overlying  a  Q  half-space  of  2COO.  By  fitting  their  data  from  2  to  28 
Hz,  they  found  a  satisfactory  fit,  but  required  a  frequency-dependent 
total  Q  due  to  the  velocity  structure  effects.  Our  data  here  cannot 
be  used  to  confirm  or  refute  these  conclusions,  as  our  frequency 
resolution  is  too  poor  at  low  frequencies.  Our  solutions  do  not 


require  frequency-dependent  Q  to  satisfy  the  data.  Further,  Lewis 
ard  Jung  [1989]  conclude  that  the  combined  effects  of  the  shallow 
highly  attenuating  layer  and  velocity  structure  limit  standard  seismic 
refraction  data  to  be  below  30  Hz.  The  modeling  presented  here 
(Figure  9)  requires  only  that  path-averaged  Q  of  the  crust  be  less 
than  about  100  to  limit  tne  frequency  content  of  typical  seismic 
refraction  data. 


Conclusions 

The  results  of  the  first  determination  of  compressional  wave 
attenuation  in  the  uppermost  oceanic  crust  revealed  that  attenuation 
can  largely  be  described  by  a  linear  frequency  dependence  of 
attenuation  and  a  frequency- independent  component  of  attenuation. 
Values  of  Q  vary  between  20  and  50  for  the  uppermost  650  m  of  the 
crust.  The  attenuation  is  significantly  higher  than  those  found  by 
modeling  data  with  synthetic  seismograms  but  is  in  general  agree¬ 
ment  with  recent  laboratory  data  [Wepfer  and  Christensen,  1987, 


I - I 

00  0.2 

Time,  seconds 

Fig.  10.  Equal  amplitude  plots  of  shots  recorded  b>  SFH8  oriented  prallcl 
and  perpendicular  to  the  ridge  crest  at  about  2  km  range.  The  shot  perpen¬ 
dicular  to  the  ridge  crest  lies  outside  the  receiving  array  and  could  not  be 
adequately  located  for  inclusion  in  the  data  inversion.  The  disparity  in 
amplitudes  suggests  an  azimuthal  variation  in  attenuation,  possibly  due  to 
vertically  oriented  cracks  and  fissures  aligned  parallel  to  the  ridge  cresL  A 
test  of  the  data  revealed  no  azimuthal  variation,  largely  due  to  the  paucity  of 
ahcts  oriented  perpendicular  to  the  ridge  crest. 
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this  issue].  Bused  upon  modeling  of  typical  seismic  refraction  data 
using  Lhe  sonar  equation,  the  results  indicate  that  our  values  may  be 
typical  to  higher  than  normal  values  of  attenuation  for  oceanic  crust. 
It  is  unlikely,  however,  that  values  of  Q  of  20  are  present  throughout 
the  whole  crust. 

Our  data  do  not  correspond  to  any  rate-scale  geological  interpreta¬ 
tion  of  the  uppermost  crust,  but  this  may  be  a  result  of  an  overly 
smooth  velocity  depth  function,  which  affects  the  attenuation 
models.  Our  data  do  indicate  a  degree  of  heterogeneity  within  the 
oceanic  crust.  The  magnitudes  ar.d  scales  of  such  heterogeneities 
cxmoi  be  determined,  bui  the  spatial  scales  are  likely  to  be  of  the 
order  of  hundreds  of  meters.  No  azimuthal  variation  of  attenuation 
was  evident,  although  the  data  set  is  somewhat  inadequate  to  con¬ 
clude  this  result. 

Further  experiments  and  results  will  be  necessary  to  determine 
whether  our  results  are  typical  or  anomalous  for  young,  upper 
oceanic  crust.  Until  such  time,  however,  a  healthy  caution  to  using 
the  values  of  Q  determined  from  synthetic  seismograms  is  war¬ 
ranted. 
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